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ABSTRACT 

We present the surface density of luminous active galactic nuclei (AGN) associated with a uniformly 
selected galaxy cluster sample identified in the 8.5 square degree Bootes field of the NOAO Deep Wide- 
Field Survey. The clusters are distributed over a large range of redshift (0 < z < 1.5) and we identify 
AGN using three different selection criteria: mid-IR color, radio luminosity, and X-ray luminosity. 
Relative to the field, we note a clear overdensity of the number of AGN within 0.5 Mpc of the cluster 
centers at z > 0.5. The amplitude of this AGN overdensity increases with redshift. Although there are 
significant differences between the AGN populations probed by each selection technique, the rise in 
cluster AGN surface density generally increases more steeply than that of field quasars. In particular, 
X-ray selected AGN are at least three times more prevalent in clusters at 1 < z < 1.5 compared to 
clusters at 0.5 < z < 1. This effect is stronger than can be explained by the evolving median richness 
of our cluster sample. We thus confirm the existence of a Butcher-Oemler type effect for AGN in 
galaxy clusters, with the number of AGN in clusters increasing with redshift. 

Subject headings: galaxies: active - galaxies: clusters: general - galaxies: evolution - infrared: galaxies 
- X-rays: galaxies - radio continuum: galaxies 



1. INTRODUCTION 

In the local universe, AGN are known to be much less 
common in rich galaxy cluster environments as compared 
to the field. Optical observations dating back more than 
three decades show that while 5% of massive field galax- 
ies show spectroscopic signatures of nuclear activity, only 
1% of the corresponding p opulation in local ga laxy clus- 
ters show such signatures ((Dressier et al. 198 5]). Assum- 
ing that the M % — a relation fe.g.. lMagorrian et al.lll998T ) 
holds in rich environments, this suggests that the super- 
massive black holes hosted by massive cluster galaxies 
are currently quiescent or their nuclear activity is opti- 
cally obscured. In addition, the intracluster plasma in 
rich galaxy clusters is maintained at high temperatures; 
AGN feedback has been suggested as the most natural 
heat source by numerou s authors (e.g.. lMcNamara et al.1 
l200CHBirzan et alj|2004l ). This suggests that the activity 
level of AGN in cluster galaxies, and their influence on 
the formation and evolution of clusters, may have been 
higher in the past. 

Early studies of AG N in clusters largely relied on op- 
tical spectroscopy (e.g.. lBahcall et aLlll969ft . a technique 
which is notoriously insensitive to both obscured and 
low luminosity AGN. In contrast, multi-wavelength ap- 
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proaches provide a significantly more thorough census of 
AGN. For example, Compton-thick, luminous AGN will 
be faint at the relatively low (< 10 keV) energies probed 
by Chandra and XMM/Newton and would be unlikely to 
show a UV-excess or spectroscopic signatures of nuclear 
activity at cosmological distances. Such sources would 
be easier to identify at radio or mid-IR energies, though 
only a ~ 10% of AGN are radio-loud. We briefly review 
selection of AGN at non-optical wavelengths, summariz- 
ing previous work using these techniques to study the 
cosmic evolution of AGN in galaxy clusters. 

X-Ray Selected AGN. In the past decade, sensitive X- 
ray observatories have provided an efficient method to 
identify AGN. Several studies conducted on individual 
clusters or small cluster samples have shown an excess of 
X-ray point source s in galaxy clusters (jCappi et alJl200it 
lEckart et ai1l2006l ). Martini et al. (2002, 2006, 2007) 
study samples of X-ray sources near low-redshift galaxy 
clusters with spectroscopic follow-up to confirm cluster 
membership. Similar to field X-ray source studies, this 
work has shown that X-ray emitting galaxies in clusters 
often do not show AGN signatures in their optical spec- 
tra, implying that a significant fraction of cluster nuclear 
activity is optically obscured. C onsidering a sample o f 
four galaxy clusters at z ~ 0.6, lEastman et al.1 (2007) 
show that the X-ray AGN are more common at z ~ 0.6 
than at z ~ 0.2, implying that the cluster AGN fraction 
increases rapidly with redshift. 

Other groups have foregone spectroscopic observations 
and have looked for enhancements in the radial distri- 
bution of X-ray sources aro und clusters. For example, 
iRuderman fc Ebelind (|2005l) der ive the surface density 
profile of X-ray sources around a sample of 51 galaxy 
clusters at 0.3 < z < 0.7 from the Chand ra MAssive 
Cluster Survey fMACS: lEbeling et alJl200l . They find 
that X-ray sources show a very pronounced central spike 
(r < 0.5 Mpc) and a second excess at ~ 2.5 Mpc relative 
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to the cluster centers. They suggest the former is from 
nuclear activity triggered by infalling galaxies approach- 
ing the central giant elliptical galaxy, while the second 
peak, which is most pronounced in virialized clusters, 
is due to galaxy merg ers at the cluster-field interface. 
iBignamini et al] ([2008) find an excess of X-ray point 
sources near the centers of optically-selected clusters at 
0.6 < z < 1.2. 

Radio Selected AGN. Over the past 50 years there have 
been numerous studies of the environments of power- 
ful radio galaxies which are now well known to often 
be associated with the centr al member of galaxy clus- 
ters, not only at low re dshifts ([Matthews et al.lll964f) but 
also at high redshifts (lDickinsonlll993t iStern et al.ll2003t 
iMilev et al.ll2004t IVenemans et al.ll2005ll2007l Galametz 
et al., in pr ep.). Correlating the Slo an Digital Sky Sur- 
vey (SDSS; iGunn & Weinberg! [l995h with the Faint Im- 
ages of the Radio Sky at Twenty cm survey (FIRST; 
iBecker et all 11995( 1 . iBest et al.l ([2005D show that there 
is a strong tendency for low-redshift radio-loud galaxies 
to be found in rich environments. Lin & Mohr (2007) 
derive the radio-loud AGN density profile from a sam- 
ple of 573 nearby (z < 0.2) X-ray selected galaxy clus- 
ters and find a c lear excess o f AGN in the center of the 
galaxy clusters. ICroft et all (|2007l ). studying brightest 
clusters galaxies (BCGs) at z < 0.3 from a sample of 
13, 240 galaxy clusters, report that 19.7% of BCGs are 
detected in the radio, with the majority of these (92.8%) 
being radio-loud (Li.4GHz > 10 23 W Hz -1 ). In con- 
trast to these studies of X-ray and radio-selected AGN in 
galaxy clusters, optically-selected, luminous quasars ex- 
hibit clustering comparable to that of local field galaxies, 
and therefo re primarily reside in lower mass da rk matter 
halos (e.g JCroom et alJl2005HCoil et al.ll2007[ l. 

Mid-IR Selected AGN. AGN can also be efficiently 
identified using mid-infrared colors (jLacy et al.l 120041 : 
IStern et all 12005( 1 . While the optical to mid-infrared 
(A < 5/im) continuum of normal galaxies is dominated 
by the composite stellar black body emission and peaks 
at approximately 1.6/im, the optical to mid- infrared 
continuum of active galaxies is dominated by a power 
law. Consequently, with sufficient wavelength baseline, 
one can easily distinguish AGN from stellar popula- 
tions. Little work has been done to date on the preva- 
lence of mid-IR selected AGN in clusters, though Hickox 
et al. (submitted) shows that low redshift (z < 0.8) 
mid-IR selected AGN have weaker clustering properties 
(e.g., smaller correlation lengths) than either X-ray or 
radio-selected AGN. At high redshift (z ~ 2), how- 
ever, various mid-infrared AGN samples show strong 
clustering, suggesti ng that they are preferentially associ- 
ated with clusters (iBrodwin et al.l l2008t iDev et"aHl2008t 
iMagliocchetti et al.ll2008D . 

Previous studies of AGN activity in galaxy clusters 
have typically focussed on a single selection technique, 
generally over a limited redshift range and with a lim- 
ited cluster sam ple size. Combining t he IRAC Shal- 
low Survey (ISS; lEisenhard t et al.1l2004t) and the Chan- 
dra X Bootes Survey ([Murray et alj|2005h . iGoriian et al.l 
(2008) compare mid-IR and X-ray AGN selection and 
find that only 28% of XBootes sources are selected by 
the IStern et al.l ([2005( 1 mid-infrared criteria. On the 
other hand, a significant number (43%) of AGN selected 
in mid-infrared have an X-ray counterpart in XBootes 



(V.Gorjian, private communication). Furthermore, only 
~ 10% of optically-selected luminous quasars are radio- 
loud (e.g.. IStern et al.ll2000l ). These results imply that 
to obtain the fullest census of nuclear activity, multiple 
AGN selection techniques are required. 

In this paper, we present the surface density of AGN 
selected by three techniques (mid-infrared color, radio 
luminosity, and X-ray luminosity) in a large sample 
of uniformly- selected galaxy cluste rs identified in the 
Bootes field (|Eisenhardt et all [2008). The next section 
(§2) describes the galaxy cluster and AGN samples. We 
describe the radial surface density profile of AGN in 
galaxy clusters in §3 and discuss the results in §4. We as- 
sume a ACDM cosmology with Hq = 70 km s _1 Mpc -1 , 
fi TO = 0.3 and fl\ = 0.7. 

2. THE GALAXY CLUSTER AND AGN SAMPLES 

The galaxy cluster sample we analyze has been uni- 
formly selected from the 8.5 square degree ISS field in 
Bootes survey, which has been observed by a wide ar- 
ray of facilities at multiple wavelengths. The following 
sections describe first the cluster sample and then AGN 
samples identified in the Bootes field at mid-infrared, ra- 
dio and X-ray energies. Note that in the later analysis 
(§3), we will use these AGN samples, generally from flux- 
limited catalogs, to study the cosmic evolution of AGN 
in galaxy clusters. To ensure uniformity in the analysis, 
this will entail creating luminosity-selected sub-samples 
(§3). 

2.1. The IRAC Shallow Survey cluster sample 

lEisenhardt et alj ([2008( 1 present a sample of 335 galaxy 
clusters identified in the Boo tes field of the NOAO Deep 
Wide-Field Survey (NDWFS;[jannuziTDe3[i99i). The 
clusters, reaching out to redshift z = 2, were identi- 
fied by combining ISS with By/RI imaging from the 
NDWFS and JK S imaging from the FLAM I NGO S 
Extragalactic S urvey (FLAMEX; lElston et all 12006( 1 . 
iBrodwin et al.l (|2006f ) calculated photometric redshifts 
for 4.5/im-selected sources in the ISS. Clusters were 
then identified using a wavelet algorithm which iden- 
tified three-dimensional overdensities of galaxies (in 2 
angular dimensions) with comparable photometric red- 
shift s. The clusters have average halo masses of ~ 10 14 
Mq ([Brodwin et al.ll2007l ). The cluster centers were de- 
termined from their peaks in the wavelet detection map. 
Since this map is based on the positions of a limited num- 
ber of galaxies per cluster, the cluster center accuracy is 
expected to be comparable to the 12" resolution of the 
wavelet map. 

An estimate of the cluster redshifts was obtained from 
the photometric redshift probability distribution of can- 
didate cluster members. At z < 0.5, 67 of the 91 clusters 
candidates are confirmed by the AGN and Galaxy Evo- 
lution Survey (AGES; Kochanek et al., in prep.), which 
provides spectroscopic redshifts for ~ 17000 objects in 
the ISS. To date, 12 of the z > 1 cl uster candidate s have 
been confir med spectroscopicallv ([Stanford et al.l 120051: 
Elston et all2006t ; IBrodwin et al.ll2~006l : lEisenhardt et all 
2008). Where available, we use the spectroscopically de- 
termined cluster redshifts. For the remaining cluster 
sample, we use the best fit photometric redshift which 
includes a small, empirically-determined offset relative 
to the mean photometric redshift of the candidate clus- 
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ter members. This large, uniformly selected galaxy clus- 
ter sample is spread over a large range of cosmic time 
and thus provides a powerful tool for studying the de- 
pendence of the AGN distribution on both redshift and 
environment. 

2.2. Mid-IR selected AGN 

Using mid-infrared colors, provided by the ISS, of 
nearly 10, 000 spec t roscop ically identified sources from 
AGES. [Stern et all (|2005l ) show that AGN can robustly 
be identified from their IRAC colors (see also Lacy et al. 
2004; Donley et al. 2008). Simple mid- infrared color cri- 
teria select over 90% of the spectroscopically identified 
type 1 AGN and 40% of the type 2 AGN. The selection 
relies on the fact that while stellar populations fade long- 
ward of their 1.6/im peak, luminous AGN exhibit red 
power-law emission throughout the mid-infrared. No- 
tably, this power-law long wavelength emission is rela- 
tively immune to gas and dust absorption, making mid- 
infrared AGN selection sensitive to both obscured and 
unobscured luminous AGN, relatively independent of 
redshift. We note however that this selection fails for 
low-luminosity AGN where the host galaxies contributes 
much of the mid-IR flux (Eckart et al., in prep.). We ex- 
tract all 2262 sources selected by the Stern et al. (2005) 
criteria that are well detected (5a) in all four IRAC 
bands. For a typical AGN spectral energy distribution, 
the 5.8/im band is the least sensitive, with a 5a limit- 
ing depth of 51/xJy (3" diameter aperture). We find a 
mean surface density of 283 mid-infrared selected AGN 
candidates per square degree in the Bootes field. 

2.3. Radio selected AGN 

Several shallow (> lmJy) radio surveys have covered 
the ISS field at 1.4 GHz, including the NRAO VLA 
Sky Survey (NVSS; Condon et al. 1998) and FIRST. 
More recently, the Westerbork Synthesis Radio Tele- 
scope (WSRT) has covered approximately 6.68 square 
degre es of Bootes, reachin g a la limit of 28/xJy at 1.4 
GHz (|de Vries et al.|[200l . The final catalog of 5a de- 
tected sources contains 3172 objects 1 , corresponding to 
a surface density of 453 radio sources per square degree. 

The brightest radio sources in the sky tend to be AGN, 
with their radio emission dominated by synchrotron ra- 
diation. However, sufficiently deep radio surveys are also 
able to detect radio emission associated with stellar pro- 
cesses such as supernova explosions and compa ct stellar 
remna nts, abundant in star-forming galaxies. iCondonl 
( 1992) shows that for these star forming galaxies, the star 
formation rate (SFR) is proportional to the radio emis- 
sion. Following Condon et al. (1992), Serjeant, Gruppi- 
oni & Oliver (2002) obtain Li. 4GHz (erg s" 1 Hz" 1 ) = 8.2 
x 10 27 x SFR(M© yr" 1 ). Submillimeter galaxies corre- 
spond to the objects with the highest known SFR, reach- 
ing SFR - lOOOAf© yr" 1 (jPope et al.ll200l . We conser- 
vatively assume that no galaxy has a SFR greater than 
3000A/© yr" 1 and thus any radio source with a 1.4 GHz 
luminosity density greater than 2.46 x 10 31 erg s _1 Hz -1 
must host an AGN (see §3). Given the depth of the 
WSRT observations, all radio-loud AGN in the Bootes 
field should be detected out to redshift z ~ 3.7. 

1 The catalog as we ll as the mosaic can be obtained through 
anonymous ftp to ftp://ftp.nfra.nl/pub/Bootes 



2.4. X-ray selected AGN 

XBootes is the largest contiguous survey with the 
Chandra X-Ray Observatory. It consists of 9 square de- 
grees of imaging of the NDWFS Bootes field, comprised 
of 126 5 ks pointings, with a corresponding depth of 
- 7.8 x 10" 15 ergs cm" 2 s" 1 in the 0.5-7 keV band. 
The X-ray catalog consists of 3293 point sources with 
four or more X-ra y counts (366 sources pe r square de- 
gree; iKenter et aLH2005t iBrand et all I2006D and is pub- 
licly available 2 . At the shallow depth of these data, the 
majority of the X-ray sources are expected to be AGN. 
Indeed, based on AGES op tical spectra for 892 XBootes 
sources (IKenter et al.ll2005h . AGN represent at least 69% 
of the sources with spectra. Only 3% of the sources 
are identified as stars. The remaining 28% lack obvious 
AGN features in their optical spectra but are suspected 
to be AGN based on their X-ray luminosities. XBootes 
is only sensitive to starbursts out to z — 0.12 (e.g., see 
iGoriian et alj [2008) . Since only one of the galaxy clus- 
ters in the ISS sample is at z < 0.12, we consider the 
possible influence of X-ray detected starbursts negligible 
in this work. 

3. RADIAL SURFACE DENSITY PROFILE OF AGN 

In order to study the cosmic evolution of AGN in 
galaxy clusters, we opt for a simple empirical approach: 
we calculate the radial distribution of AGN around 
galaxy clusters as a function of cluster redshift, with 
AGN identified using the three selection criteria de- 
scribed above. We separate the ISS clusters into three 
redshift bins: < z < 0.5, 0.5 < z < 1 and 1 < z < 1.5, 
resulting in 91, 140 and 79 clusters per bin, respectively, 
for the mid-IR and X-Ray techniques. The remaining 25 
clusters are at z > 1.5, and are thus ignored in this work. 
The area covered by the WSRT radio data is smaller than 
the ISS and thus only 285 clusters have deep coverage at 
radio wavelengths. Of these 285 clusters, 77, 121 and 69 
are found at < z < 0.5, 0.5 < z < 1 and 1 < z < 1.5, 
respectively. 

In order to ensure that similar populations are plotted 
in each cluster redshift bin, we apply uniform luminos- 
ity cuts to the data. In detail, since most of the sources 
considered here lack spectroscopic redshifts, for a given 
cluster at redshift z, we only consider sources brighter 
than some evolving flux density such that their lumi- 
nosity would be equal to or above our luminosity limit 
assuming that the source were in the cluster. For the 
very sensitive radio survey, we apply the luminosity cut 
Si.4GHz > 2.5 x 10 40 ergs s" 1 (see §2.3) to isolate AGN 
from star-forming galaxies. For the mid-IR and X-ray 
selected AGN, we apply a luminosity cut corresponding 
to the flux limit of the corresponding surveys at z = 1 
(see §3.1 and 3.3 for details). This redshift is selected 
as a compromise between obtaining sufficient numbers of 
AGN, large numbers of clusters, and sufficient leverage 
on the cosmic time probed. Table 1 presents the num- 
ber of sources that remain for each selection technique 
and redshift bin. These numbers are based on cuts (de- 
scribed in detail in the following subsections) correspond- 
ing to the median redshift of the clusters in each bin, e.g., 
(z) = 0.38 (0.39), (z) = 0.74 (0.74) and (z) = 1.19 (1.18) 
for mid-IR and X-ray (radio) selections. Also tabulated 

2 http:/ /heasarc. nasa.gov/W3Browse/all/xbootes. html 
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are the field AGN surface densities for each technique and 
rcdshift bin, where we have simply divided the number 
of sources by the corresponding survey area. These field 
surface densities are shown as dashed horizontal lines in 
Fig. [U Finally, we note that for the X-ray and mid-IR 
selected AGN in the highest redshift bin (1 < z < 1.5), 
this approach only provides a lower limit to the surface 
density of AGN down to the luminosity limit considered 
in the lower redshift bins. 

The surface density distributions for AGN in the 
galaxy clusters are calculated as follows. We calculate 
the angular separation from each ISS cluster center to 
each mid-infrared, radio and X-ray selected AGN can- 
didate, subject to the flux cuts discussed above. AGN 
candidates are identified out to 5' from the cluster cen- 
ters, and we bin the radial distances to compute the sur- 
face density profiles of AGN (0.5' bins), presented in 
Fig. [TJ As expected, at large radii the derived surface 
densities asymptote to the field AGN surface densities. 
Throughout this paper, we adopt Poissonian errors for 
AGN counts and use the iGehrelsl (|1986h small numbers 
approximation for Poisson distributions. Fig. [T] presents 
the la errors for both the AGN counts and the field sur- 
face densities for each technique and redshift bin. 

As a check on this method, we also calculated the AGN 
surface density distribution at random locations in the 
field. We created 40 mock cluster samples by simply 
offsetting all of the cluster positions by medium-scale 
offsets — larger than the typical cluster size, but small 
enough that few clusters would fall outside of the survey 
region. Deriving the radial surface density profile around 
the mock cluster samples, we consistently find no signif- 
icant variation relative to the field surface densities. 

In contrast to the mock sample, an excess of AGN is 
found near the centers of clusters, most prominently at 
z > 0.5. The next sections detail the results obtained 
for each of the three AGN selection criteria considered. 
We note that we have verified that spatially proximate 
clusters are not biasing our results. Since the overdensi- 
ties described below are only found at small radii (< 1 
Mpc), the number of AGN counted as potential central 
members to more than one cluster is negligible (null for 
r < 0.5 Mpc, and < 1% for 0.5 < r < 1 Mpc). 

3.1. Density profile of mid-IR selected AGN 

The radial surface density profile for mid-IR selected 
AGN is shown in Fig. [1] (top row), separated into the 
three redshift bins described above. The 5.8/im limit- 
ing flux density of the ISS is 51/iJy (5a). Assuming a 
pure power law spectrum for the AGN, the correspond- 
ing luminosity density is L v — A-ko1 2 l S u / (1 + z)' 1- "' 
where dj, is the luminosity distance and a is the spec- 
tral index (S v oc v a ). This assumption of a power- 
law is born out by both the mid-IR color-color plots of 
broad-lined AGN used to develop the mid-IR AGN se- 
lection criteria (e.g., Lacy et al. 2004, Stern et al. 2005) 
and bro ad-band composite A GN spectral energy distri- 
butions (|Richards et alj 120061 ). We adopt a mean spec- 
tral index of 0.73, as given by iStern et ail (|2005l ). The 
depth of the ISS 5.8/im data corresponds to a 5.8/im 
luminosity of 2.2 x 10 44 ergs s _1 at z — 1, which we 
adopt as our luminosity cut. Thus for any given cluster 
we only consider mid-IR selected AGN whose luminosity 
(and corresponding flux density) is above this value un- 



der the assumption that they are at the cluster redshift. 

The z < 0.5 galaxy clusters do not show any signif- 
icant enhancement of mid-IR selected AGN relative to 
the field. In the 0.5 < z < 1 bin, however, a weak over- 
density of mid-IR selected AGN is observed at the cluster 
center (r < 0.3 Mpc). Throughout the paper, the signif- 
icance of overdensities have been derived as follows. We 
first subtract the number of AGN in a given aperture by 
the field AGN population to a given (redshift-dependent) 
flux density limit. This field-corrected number of AGN 
is then referen ced to the expected field counts using the 
IGehrelsl (|T986) small numbers approximation for Pois- 
son distributions. The central excess of 0.5 < z < 1 
mid-IR selected AGN is thus significant at the 2a level. 
A small excess of AGN is suggested in the center of the 
z > 1 galaxy clusters, but is not significant relative to the 
field density. Recall that since the highest redshift bin 
considers more luminous AGN than the z < 1 bins, the 
resulting, weak overdensity is a lower limit to the number 
of AGN per cluster to the luminosity limit of the z < 1 
clusters. 

3.2. Density profile of radio selected AGN 

The radial surface density profile for radio-selected 
AGN is shown in Fig. Q] (middle row), separated into 
three redshift bins. For each redshift bin, we apply the 
conservative uniform luminosity cut described in §2.3 
to remove star-forming galaxies, adopting a mean radio 
spectral index of —0.75, as given in de Vries et al. (2002). 

In the z < 0.5 bin, a weak overdensity (1.5a) of radio 
sources is observed at the cluster center. The 0.5 < z < 
1 galaxy clusters show a very clear excess of AGN in 
the center (r < 0.3 Mpc), with no excess seen at larger 
radii. The central excess is significant at the 3a level. 
The higher redshift galaxy clusters (1 < z < 1.5) also 
show a significant excess (2a) at the cluster center, albeit 
weaker and on a larger scale (r < 0.5 Mpc) relative to 
the moderate redshift galaxy clusters. 

3.3. Density profile of X-ray selected AGN 

The radial surface density profile for X-ray selected 
AGN is shown in Fig. [1] (bottom row). The limiting 
flux of XBootes in the 0.5 — 7.0 keV range is ~ 7.8 x 
10~ 15 ergs cm~ 2 s^ 1 , where the XBootes flux values have 
been calculated from the counts assuming an X-ray AGN 
power-law spectrum with photon index T = 1.7 (Kenter 
et al. 2005). While the shallow depth of the XBootes 
survey implie s that any source de tected at z > 0.12 must 
be an AGN (|Goriian et al.|[2008| ). we only plot sources 
whose X-ray luminosity exceeds the limiting sensitivity 
of XBootes at redshift unity. The same photon index 
value is adopted to calculate the X-ray luminosities. The 
luminosity cut applied for X-ray selection is therefore 3.3 
x 10 43 ergs s _1 in the 0.5 — 7.0 keV range. As emphasized 
for the mid-IR selected AGN, the lower, right panel of 
Fig. [1] is a lower limit to the surface density of X-ray 
selected AGN at the luminosities plotted in the first two 
panels. 

No significant enhancement of the X-ray selected AGN 
is observed for the lowest redshift (z < 0.5) galaxy 
clusters. A small central overdensity is suggested for 
0.5 < z < 1 galaxy clusters (1.2tr level; r < 0.5 Mpc). On 
the other hand, the 1 < z < 1.5 galaxy clusters show a 
modest excess (2a) at the cluster center (r < 0.25 Mpc). 
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We note that, in principle, X-ray point sources are 
harder to identify in galaxy clusters with significant dif- 
fuse X-ray emission. The intergalactic medium in such 
clusters would tend to decrease the apparent AGN sur- 
face density for lower redshift clusters in a shallow survey 
such as XBootes. X Bootes identifies 43 extended sources 
(|Kenter et al.ll2005l ) of which 38 are covered by ISS. The 
size of these extended sources, derived from Gaussian 
fits to their source profiles, is (except for one source) less 
than 0.5'. The overdensity of X-ray point sources that 
we observe in our galaxy cluster sample is clearly de- 
tected within 0.5 Mpc (r ~ 1.5' at z = 0.5) of the cluster 
center. The brightness of an unresolved cluster mem- 
ber with an X-ray luminosity equivalent to our thresh- 
old would, on average, be five times brighter than that 
of the extended sources. Only in four clusters would a 
threshold source have a flux comparable to the extended 
source (e.g., brightness contrast < 2). Given these small 
numbers, we are confident that extended, diffuse X-ray 
emission from rich clusters is not strongly biasing our 
results. 

4. DISCUSSION 

This paper presents the first comprehensive study of 
the radial density profile of AGN for a uniformly-selected 
cluster sample distributed over a large fraction of cosmic 
time. Only the radio-selected AGN show a weak (1.5a) 
overdensity in the lowest redshift (z < 0.5) bin. An ex- 
cess of AGN is observed near the center (r < 0.5 Mpc) 
of galaxy clusters at z > 0.5, and all three AGN selec- 
tion criteria show a pronounced enhancement in AGN 
activity at 0.5 < z < 1 relative to z < 0.5. The excess 
is even observed for the highest redshift galaxy clusters 
(1 < z < 1.5), though we are only plotting to the flux 
limits of the input mid-IR and X-ray catalogs at these 
high redshifts, implying a less sensitive luminosity limit 
than in the z < 1 clusters. We therefore show that AGN 
activity in galaxy clusters increases with redshift. 

We first compare the luminosity cuts applied for each 
selection criteria to see how the AGN samples relate 
to one another. The luminosity cut applied for the 
mid-infrared sele ction is 2.2 x 10 44 ergs s _1 at 5.8/im. 
lElvis et all (|1994D provides a mean spectral energy distri- 
bution (SED) for quasars derived from broad-band pho- 
tometry covering the full electromagnetic spectrum from 
radio to X-rays. Using this composite SED for radio- 
loud quasars, the luminosity cut at 5.8/im corresponds 
to a luminosity cut of ~ 2.2 x 10 43 ergs s _1 in the 0.5 — 7 
keV X-ray band and 2.2 x 10 41 ergs s" 1 at 1.4 GHz. Our 
X-ray luminosity cut is 3.3 x 10 43 ergs s _1 which is com- 
parable to the depth of the mid-infrared selection. A 
more recent d etermination of the m ean quasar SED was 
performed by iRichards et al.l (|2006l ) and gives an X-ray 
luminosity five times smaller at a given mid-infrared lu- 
minosity compared to the Elvis et al. (1994) template. 
This would imply that our X-ray selection only identifies 
AGN five times more powerful than the mid-infrared se- 
lection. Our radio luminosity cut is 2.5 x 10 40 ergs s _1 
at 1.4 GHz. Therefore, the AGN radio selection depth is 
reaching almost ten times the depth of the mid-infrared 
selection, and ten to fifty times the depth of the X-ray se- 
lection, depending on which quasar composite template 
is adopted. This deeper radio sensitivity is a plausible 
explanation for why only the radio-selected AGN show 



an overdensity in the lowest redshift galaxy clusters. 

One concern is that due to the larger volumes acces- 
sible at higher redshifts, our sample of clusters becomes 
progressively richer with redshift. In this case, the ris- 
ing number of AGN seen in the higher redshift clusters 
could simply trace the rising number of massive cluster 
members in the higher redshift clusters, with a constant, 
non-evolving fraction of massive galaxies being active. 
In order to study this effect, we use the photometric red- 
shifts derived from the 4.5/xm ISS sample (see §2.1) to 
determine the mean number of cluster members brighter 
than 0.5L* as a function of cluster redshift. We con- 
sider the observed 3.6/im flux density, converted to 1.6/xm 
rest- frame luminosity assuming azy = 3 formation red- 
shift and passive evolution (e.g., see lEisenhardt et al] 
2008). We only consider members brighter than 0.5L* 
within the inner 0.5 Mpc of the clusters. As expected, 
we do indeed find clusters at z ~ 0.3 are, on average, 
slightly less rich than clusters at z ~ 0.8. Clusters at 
< z < 0.5 have an average of 6.5 members brighter 
than 0.5L* while 10.0 members are found, on average, in 
clusters at 0.5 < z < 1. At redshifts above 0.8, the num- 
ber of cluster members in the inner 0.5 Mpc and brighter 
than 0.5L* actually begins to decrease, with 5.7 galaxies, 
on average, for clusters at z > 1. This is not a sensitivity 
issue since the 90 seconds of combined exposure of ISS 
provides suffic ient sensitivity to detec t evolving L* galax- 
ies to z — 2 (|Eisenhardt et al.| [2008). This decrease in 
cluster member numbers instead could reflect that clus- 
ters at z > 0.8 are still in the process of collapsing, that 
the clusters members are merging, or some combination 
thereof. However, the variation with redshift is only at 
the 50% level, while the cosmic evolution of AGN ac- 
tivity is a factor of several more significant (with some 
variation depending on which AGN selection criterion 
has been used). For the following discussion, we there- 
fore omit the modest systematic dependence of cluster 
richness on redshift inherent to the Bootes cluster sam- 
ple. We do note, however, that this trend would serve to 
make the fraction of active galaxies in the central regions 
of z > 1 clusters even higher relative to the fraction at 
z < 0.8, since the number of AGN is rising while the 
number of luminous cluster members is falling. 

It is well known that AGN are more common at high 
redshift. Indeed, there is a strong increase in the quasar 
density with redshift from the local unive rse out to z ~ 2 , 
a result that dates back four decades (|Schmidtl fl968). 
Since then, numerous studies have determined the quasar 
luminosity function (QLF) for radio (e.g., Dun lop & Pea- 
cock 1990), optical (e.g ., ICroom et al.1 12004). infrared 
(e.g.. iBrown et al"]l2006l) and X-ray AGN samples (e.g., 
lUeda et al.l I2003Q . Eastman et all (|2007t ) find that the 
cluster AGN fraction at z ~ 0.6 is approximately a factor 
of 20 greater than the cluster AGN fraction at z ~ 0.2. 
In comparison, this increase is only 1.5 (3.3) for the field 
sample of AGN more luminous than L x > 10 42 (10 43 ) 
ergs s _1 over this redshift range. In order to compare the 
AGN overdensities associated with galaxy clusters, seen 
in Fig. [U to the field AGN population, we derive aver- 
age cluster AGN volume densities as follows. We count 
the number of sources projected within 0.25 Mpc of the 
cluster centers for the various redshift bins and selection 
criteria, and then field-correct this number by subtract- 
ing the expected number of sources over this same area 
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based on the field surface density. We assume any iden- 
tified overdensity is associated with the cluster (e.g., is 
neither foreground nor background). For the sake of sim- 
plicity, we assume a spherical distribution at the cluster 
rcdshift to derive a volume density from the measured 
(field-corrected) surface density. Fig. [2] presents this 
field-corrected volume density of AGN within 0.25 Mpc 
of the galaxy clusters in the three redshift intervals and 
for the three AGN selection criteria. We also plot the 



Croom et alJ ( 


2004) 


Croom et al.l ( 


2004) 



shift range using the QLF derived by 
from the 2dFQZ project 3 . At z = 1 

find that L* AGN = 1.6 x 10 45 ergss" 1 , corr esponding to 
M B = -24.5. Using the lElvis etaTI (|1994D quasar tem- 
plate, the mid-IR and X-ray depth of our survey corre- 
spond to optical luminosities of ^agn + I {Mb = —23.5). 
Therefore, we plot in Fig. [5] the cosmic evolution of 
optically-selected field quasars brighter than £ AGN + 1 
(A1 B = -23.5) and L\ GN (M B = -24.5). We find a 
more pronounced cosmic evolution in the relative den- 
sity of cluster AGN compared to the field from z = 
to 1.5. In particular, the field-corrected X-ray selected 
AGN density in clusters evolves at least three times more 
rapidly than that of the field population. 

This work shows that a higher surface density of AGN 
is observed in galaxy cluster centers than in the field and 
that this AGN overdensity increases with redshift out to 
z > 1. We therefore confirm a Butc her-Oemler type ef- 
fect for AGN in galaxy clusters (e.g. JMartini et al.ll2007l) 
as the number of AGN in clusters increases with redshift 
similar to the way that higher redshift galaxy clusters 
have larger numbers of blue, star-forming galaxies. 

We find that up to ten percent of the galaxy clusters 
have at least one AGN near their center (r < 0.25 Mpc). 
Table 2 summarizes the percentage of galaxy clusters 
with an AGN for each selection technique and redshift 
bin. These fractions are field-corrected. Considering the 
luminosity cuts applied for each selection, only 3% of 
galaxy clusters at z < 0.5 have a radio-selected AGN 
near their center (r < 0.25 Mpc) and neither mid-IR nor 
X-ray selected AGN are found. We find that the fraction 
of galaxy clusters hosting such an AGN evolves with red- 
shift. By z ~ 1.5, up to 10% of our clusters sample host 
at least one X-ray emitting AGN near their center (9% 
for radio-selected AGN). 

Finally, we note that the current analysis is rather 
simple, using only the apparent two dimensional distri- 
bution of sources across the sky to address the evolv- 
ing role of AGN in clusters. In principle, a three di- 
mensional analysis, accounting for source redshifts (ei- 
ther spectroscopic or photometric), would significantly 
lower the field counts in each redshift bin and provide a 
more robust measurement of the AGN Butcher-Oemler 
effect. We have opted for the former analysis as it re- 
lies only on the simple observables of position and flux, 
and is thus model free. Brodwin et al. (2006) derive 
photometric redshifts for the Bootes field AGN with an 
accuracy a = 0.12(1 + z), with a 5% outlier fraction. 
This is approximately half the precision that can be de- 

3 Croom et al. (2004) present various forms of the empirical 
QLF. We use the version which fits the QLF with a second-order 
polynomial to describe the luminosity evolution (Table 5 of that 
paper). 



rived for normal galaxies. We are currently obtaining 
spectroscopic observations of candidate active members 
of clusters. Future work will use these results to bet- 
ter constrain the pronounced cosmic evolution of AGN 
activity in rich environments. 

5. SUMMARY 

The spatial distribution of AGN has been derived for 
a well-defined sample of galaxy clusters at < z < 1.5. 
The AGN are selected using three different selection tech- 
niques (mid-IR color, radio luminosity, and X-ray lumi- 
nosity) that isolate distinct types of active galaxies to 
ensure the fullest selection of AGN in our galaxy cluster 
sample. An overdensity of AGN is found near the cen- 
ter (r < 0.5 Mpc) of galaxy clusters at z > 0.5 and this 
AGN excess increases with redshift. We confirm that the 
rising number of AGN is not simply due to the clusters 
becoming richer with redshift with a non-evolving frac- 
tion of cluster members being active. It is well known 
that AGN are also more common at high redshift. We 
therefore also confirm that the observed increase in red- 
shift is more pronounced in the galaxy cluster sample 
than in the field. We thus find a Butcher-Oemler type 
effect for AGN in galaxy clusters. 
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TABLE 1 

Number and surface density of Bootes AGN candidates after luminosity cuts. 



Selection 


Luminosity limit 
(ergs s _1 ) 


Area 
(deg 2 ) 


Total number 
N (E, arcmin -2 ) 


< z < 0.5 
Nagn (Sagn) 


0.5 < z < 1 
Nagn (^agn) 


1 < z < 1.5 
Nagn (^agn) 


Mid-IR 

Radio 

X-ray 


5 5 .8 M m > 2.2 X 10 44 
51.4GHz > 2.5 X 10 40 
Sb.S-TkeV > 3-3 X 10 43 


8.0 
6.68 
9.0 


2262 (0.079) 
3172 (0.132) 
3293 (0.102) 


177 (0.006) 
229 (0.009) 
205 (0.006) 


1570 (0.054) 
591 (0.025) 
2102 (0.065) 


2262 (0.079) 
1394 (0.058) 
3293 (0.102) 



TABLE 2 

Fraction of galaxy clusters with at least one AGN in the inner 0.25 Mpc 



Selection < z < 0.5 0.5 < z < 1 1 < z < 1.5 



Mid-IR < 1% 5 ± 2% 4 ± 3% 

Radio 3 ±2% 8 ± 3% 9 ± 4% 

X-ray < 1% 3 ± 2% 10 ± 4% 
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Fig. 1. — Radial density profile of mid-IR selected AGN (top row), radio-selected AGN (middle row) and X-ray selected AGN (bottom 
row). Galaxy clusters are divided into three redshift bins: z < 0.5 (left column), 0.5 < z < 1 (middle column) and 1 < z < 1.5 (right 
column). The global field density is shown by the red dashed line. The red dotted lines indicate the lcr variations in the field densities and 
the error bars indicate the lcr errors on the AGN counts; both assume Poissonian errors and Gehrels (1986) small numbers approximation. 
The 1 Mpc and 2 Mpc projected radii are indicated by vertical dot-dashed lines for the median redshift of the clusters plotted in each 
bin. The number of galaxy clusters used in the profile calculation is also indicated at the top right of each panel. An overdensity of AGN 
is observed in the centers of galaxy clusters at z > 0.5 while little or no variation relative to the field density is found at lower redshifts 
(z < 0.5). 
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Fig. 2. — Field-corrected volume density of AGN in galaxy clusters and comparison of its evolution relative to an optical quasar luminosity 
function. Points show the field-corrected overdensity of AGN within 0.25 Mpc of the ISS cluster centers for the three AGN criteria: mid- 
infrared (black points), radio (green triangles) and X-ray (blue squares). No mid-IR or X-ray selected AGN are found within 0.25 Mpc 
for the galaxy clusters at z < 0.5. We plot the la uncertainty of the field density as an upper limit. Lines show the relative evolution of 
quasars brighter than £^GN (Mb = —24.5, dotted red line) and L"^ GN + 1 (Mb = —23.5, dashed red line) for the field, where i^GN ' s 
calculated at z = 1. These models, which show a less dramatic evolution of AGN in the field relative to clusters, are calculated from the 
Croom et al. (2004) 2dFQZ QLF (see text for details). 



